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[57] ABSTRACT 

The sensitivity model of a complex optical system ob- 
tained by linear ray tracing is used to compute a control 
gain matrix by imposing the mathematical condition for 
minimizing the total wavefront error at the optical sys- 
tem's exit pupil. The most recent deformations or error 
states of the controlled segments or optical surfaces of 
the system are then assembled as an error vector, and 
the error vector is transformed by the control gain 
matrix to produce the exact control variables which 
will minimize the total wavefront error at the exit pupil 
of the optical system. These exact control variables are 
then applied to the actuators controlling the various 
optical surfaces in the system, causing the immediate 
reduction in total wavefront error observed at the exit 
pupil of the optical system. 

12 Claims, 4 Drawing Sheets 
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FEEDBACK CONTROLLED OPTICS WITH 
WAVEFRONT COMPENSATION 

ORIGIN OF THE INVENTION 5 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the Contractor has elected not to retain 
title. 

BACKGROUND OF THE INVENTION 

1 . Technical Field 

The invention relates to precision optical devices 15 
such as astronomical reflecting telescopes with many 
reflecting or refracting surfaces whose deformations are 
to be compensated in order to maximize image quality. 

In particular, the invention relates to a feedback control 
system which computes in a single linear computation 20 
the exact set of actuator commands which minimize the 
total wavefront error at the telescope objective (output) 
for a given set of measured deformations of the reflect- 
ing or refracting surfaces of the optical device. 

2. Background Art 25 

Image quality in powerful optical devices such as 

interstellar reflecting telescopes is impaired by deforma- 
tion of the reflecting and refracting optical surfaces in 
the telescope. Deformation of the optical surfaces due, 
for example, to temperature variations, vibration, etc., is 30 
unavoidable, and must be compensated if image quality 
is to be maximized. Such compensation ideally should 
be performed quickly in real time to permit rapid repeti- 
tive compensation in order to maintain maximum image 
quality at all times through varying conditions. 35 

Alignment of a complex optical system is a related 
problem. U.S. Pat. Nos. 4,471,447 and 4,471,448 to Wil- 
liams et al. and Williams, respectively, disclose a 
method for aligning a complex optical system in which 
the desired adjustment of all of the controlled optical 
surfaces is computed using a numerical regression esti- 
mation technique. In essence, this technique perturbs 
the alignment parameters in a model of the complex 
optical system in successive trial-and-error attempts to 45 
degrade the model and duplicate the distortion actually 
observed in the (objective) output image of the optical 
system. After many such attempts, the observed distor- 
tion is ultimately duplicated in the model, and the values 
to which the alignment parameters of the model have 
been last perturbed are used to calculate the alignment 
errors in the controlled optical surfaces of the complex 
optical systems. The alignment errors are then cor- 
rected. 

The problem with such a technique is that the numer- 55 
ical regression process is so time-consuming that its use 
for anything other than initial optical alignment is im- 
practical. During operation, the corrections would be 
performed so slowly using such a technique that they 
would not keep up with the continuing perturbations of 60 
the optical surfaces and in essence would be of no prac- 
tical effect on image quality. Thus, it has not seemed 
possible to compensate for deformation of the optical 
surfaces in a complex optical system in real time. 

A related problem is that the regressive estimation 65 
approach of such a technique by definition does not 
yield the exact answer after a limited number of itera- 
tions, and therefore provides inexact compensation for 
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the deformation or movement of the optical surfaces in 
a complex optical system. 

Accordingly, it is an object of the invention to pro- 
vide a feedback control system for maintaining optical 
alignment in a complex optical system in real time. 

It is a further object of the invention to perform a 
single rapid computation of the exact alignment correc- 
tions for all controlled optical elements in the optical 
system to minimize the total wavefront error at the 
objective (output) of the optical system. 

DISCLOSURE OF THE INVENTION 

A model of the optical system is formulated from the 
design of the system optical geometry and characteris- 
tics. This model defines the system’s objective function 
(image output) in terms of deformations (error states) of 
the various segments or optical surfaces in the system 
and the actions or control variables of all of the actua- 
tors controlling the optical surfaces. (In the preferred 
embodiment, the objective function is the total wave- 
front error squared at the exit pupil of the optical sys- 
tem, the error states are 6 degree of freedom (DOF) 
rotational and translational motions of each of the opti- 
cal surfaces and the control variables are the 3 DOF 
(tilt, tip and piston) of the actuators controlling each of 
the optical surfaces. The wavefront error for each ray is 
its optical path differene with respect to an ideal or 
unperturbed state of the system.) A sensitivity model is 
formed by taking the derivatives of the objective func- 
tion with respect to the error states and control vari- 
ables. The sensitivity model is a linear model of the 
influence of the error states and control variables on the 
objective function. In one embodiment of the invention, 
the sensitivity model is then used to compute an actua- 
tor control gain matrix by imposing the mathematical 
condition for minimizing the total wavefront error at 
the exit pupil. The most recent deformations or error 
states of the controlled segments or optical surfaces of 
the system are then assembled as an error vector, and 
the error vector is transformed by the control gain 
matrix to produce the exact control variables which 
will minimize the total wavefront error at the exit pupil 
of the optical system. These exact control variables are 
then applied to the actuators controlling the various 
optical surfaces in the system, causing the immediate 
reduction in total wavefront error observed at the exit 
pupil of the optical system. In one embodiment of the 
invention, the sensitivity matrix is updated based upon 
the changes attributable to the most recent change in 
the control variables, and the entire control process is 
then repeated. 

The advantage of the invention is that each time the 
sensors sample the deformations of all the controlled 
optical surfaces of the system, the measured error states 
are multiplied by the gain matrix in a single calculation 
to provide, in a nearly instantaneous fashion, the exact 
values of the control variable which will minmize the 
total wavefront error at the exit pupil. Thus, the inven- 
tion provides an exact control solution in nearly instan- 
taneous fashion. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Preferred embodiments of the invention are de- 
scribed in detail below with reference to the accompa- 
nying drawings, of which: 

FIG. 1 is a simplified block diagram illustrating the 
concept of the feedback control system of the invention; 
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FIGS. 2 and 3 illustrate a particular optical system in 
which the invention is employed; 

FIG. 4 illustrates the coordinate system for each 
individual optical segment employed in the system of 
FIGS. 2 and 3; 

FIG. 5 illustrates the two types of sensors employed 
in the system of FIGS. 2 and 3 to measure movement or 
deformation of each optical segment; 

FIG. 6 illustates the exact ray-trace optical geometry 
employed in formulating the sensitivity matrix in one 
embodiment of the invention; 

FIG. 7 illustrates the geometry corresponding to 
FIG. 6 employed to define perturbations in the optical 
paths in the system due to deformations of the various 
optical segments; and 

FIG. 8 illustrates the geometry of the deformation 
sensors between two adjacent segments employed in 
defining a measurement matrix representing the sensor 
response employed in accordance with one embodiment 
of the invention. 

MODES FOR CARRYING OUT THE 
INVENTION 

The broad concept of the invention is illustrated in 
FIG. 1. A complex optical system has a set of optical 
surfaces 10<z, 10 b, 10c, the optical surface 10c consisting 
of two segments 10c-l and 10c-2 and an exit pupil 12. 
Each optical surface 10 has its own set of deformation/- 
movement sensors 14 and its own set of actuators 16. 
For example, the optical surface 10 a has a set of sensors 
14 a which detect motion in six degrees of freedom by 
the optical surface 10a, namely translation along three 
orthogonal axes and rotation about each of the axes. 
The remaining sensors 14 b, 14c-l and 14c-2 perform 
corresponding functions in relation to the optical sur- 
faces 106, 10c-l and 10c-2, respectively. In a preferred 
embodiment, the optical surface 10a is controlled by a 
set of actuators 16a which control rotation about two of 
the orthogonal axes most nearly parallel to the tangen- 
tial plane of the optical surface 10 and translation along 
the remaining orthogonal axis, which is normal to the 
tangential plane. The remaining actuators 1 66, 16c- 1 and 
16c-2 perform corresponding functions in relation to the 
optical surfaces 106, 10c-l and lOc-2, respectively. 

A converter/estimator 18 converts the analog out- 
puts of the sensors 14 to digital representations of the 
actual movements by the respective optical surfaces, 
comprising deflection or deformation error data. The 
error data for the i th one of the optical surfaces 10 com- 
prises the rotations 0*0), 0y( i), B^i) about orthogonal x, 
y and z axes and the translations 8*(i), 8/i), 8/i) along 
the same axes, respectively. The error data for all opti- 
cal surfaces 10 is assembled into a single error vector x. 

A sensitivity matrix generating computer 20 retrieves 
from a memory 22 the optical geometry design data 
specifying the initial shape, optical characteristics, cor- 
rect or unperturbed location and correct orientation of 
all of the optical surfaces 10. (The data stored in the 
memory 22 may be obtained using well-known initial- 
ization techniques with a pixelated sensor at the objec- 
tive pupil 12.) The computer 20 uses this data to trace 
the path of each ray in a bundle of rays 24 received at 
the input aperture 26 of the optical system comprising 
the optical surfaces 10. From this optical ray tracing, 
the computer 20 computes the sensitivity of the optical 
path length of each ra y to each deformation or element 
in the error vector x representing all of the optical 
surfaces 10. The process by which the computer 10 
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accomplishes this is discussed in further detail in the 
Working Example given below in this specification. 
The computer forms an ordered array of all such sensi- 
tivities to produce a sensitivity or C matrix, 

5 An actuator response matrix computer 28 retrieves 
from a memory 30 the actuator design data specifying 
the location, orientation and operating characteristics of 
all of the actuators 16 controlling the optical surfaces 
10. The computer 28 uses this data to compute the exact 
10 responsivity of individual components of the error vec- 
tor x to control vectors governing individual ones of 
actuators 16. The manner in which the computer 28 
accomplishes this is described in detail below. The com- 
puter 28 forms an ordered array of these responsivities 
15 to produce an actuator response or A matrix. 

A wavefront error computer 32 computes from the C 
matrix and the A matrix certain quantities (namely, 
[A^C^CA]- 1 and A^C^C) from which the total wave- 
front error at the exit pupil 12 can be computed from 
20 the error vector. These quantities together comprise an 
actuator gain matrix and are derived by satisfying the 
exact condition for minimizing the total wavefront 
error at the objective pupil 12, as will be described 
below. An actuator command computer 34 combines 
25 these quantities (the actuator gain matrix) with the error 
vector x, producing an actuator command vector u 
which defines the simultaneous optimum actuator com- 
mands for all the actuators 16 which minimize the total 
wavefront error at the objective pupil 12. The comput- 
30 ers 32 and 34 may be thought of as a single real time 
feedback control computer 36. In fact, all of the com- 
puters 20, 28 and 36 may be implemented in a single 
computer or processor system. The individual control 
commands contained in the vector u are applied to 
35 corresponding ones of the actuators 16. Digital-to- 
analog converters (not shown) may be employed to 
transform the digital outputs of the computer 34 to the 
appropriate analog voltages required to effect the me- 
chanical control excersized by each actuator 16. 

40 The advantages of the feedback control system of 
FIG. 1 include extremely high speed: each time the 
sensors 14 sample the error states of the optical surfaces 
10, a single one-step calculation of the optimum actua- 
tor commands is performed. This calculation consists in 
45 multiplying the measured error vector by the gain ma- 
trix. Also, the result is exact: the gain matrix produced 
by the computer 32 produces an exact minimization of 
the total wavefront error at the objective pupil 12. The 
invention thus provides the first high-speed exact feed- 
50 back control for complex optical systems. 

In one embodiment of the invention, the converter- 
/estimator 18 may provide estimates of the measure- 
ments of the errors in x by the sensors 10 based upon a 
measurement matrix H whose construction and use by 
55 the converter/estimator 18 is described below in detail. 

In some implementations of the invention, the com- 
manded movement of the actuators 16 may so change 
the geometrical relations among the optical surfaces 10 
as to require that the sensitivity matrix be recomputed 
60 by the computer 20 based upon a modification of the 
initial geometry data describing the optical surfaces 10. 
The modification to the initial geometry data is per- 
formed by an optical geometry change computer 38. 
The computer 38 receives the error data from the con- 
65 verter/estimator 18 and the initial geometry data from 
the memory 22 and updates the initial geometry data 
using the error data to generate a new set of geometry 
data. The computer 38 transmits this new set of geome- 
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try data to the sensitivity matrix generating computer 
20, which generates an updated C matrix therefrom. 

WORKING EXAMPLE 

INTRODUCTION TO WORKING EXAMPLE 5 

The NASA Precision Segmented Reflector (PSR) 
program is developing technologies to enable the flight 
of large spacebome astronomical telescopes starting 
sometime around the turn of the century. Much of the 
current activity centers on the PSR Focus Moderate 10 
Mission Telescope (FMMT), a medium-sized (3.65 m 
aperture) space-based submillimeter-wave telescope 
which exemplifies many of the technical difficulties that 
must be overcome. 

The example addresses the problem of figure mainte- ^ 
nance control for the FMMT primary mirror, which is 
formed of seven separate segments. Figure maintenance 
control is used when the ideal configuration of the pan- 
els with respect to each other is known. It uses sensors 
that directly sense motions of the panels to determine 20 
their errors relative to the ideal state. These errors are 
corrected by actuators that rotate and translate the 
panels. Figure initialization control is the process of 
determining the ideal configuration for the panels and 
moving them there, using a science detector. By look- 2 5 
ing at known bright calibration stars, the panels are 
moved to maximize the image quality directly. 

Figure maintenance is used between initializations to 
hold the figure of the mirror steady while taking science 
data. It is intended to counter the long-term thermal 30 
drifts that the structure experiences. For the FMMT, 
these errors are expected to be 10-200 microns at the 
segment centers over the course of a single orbit. They 
occur on a relatively long time scale, typically tens of 
minutes. Consequently, the FMMT figure maintenance 35 
controller has a low bandwidth on the order of 0.1 Hz. 

The FMMT also experiences vibrational disturbances 
due to the presence on the structure of vibrating ma- 
chinery, such as reaction wheels and a secondary mirror 
chopping mechanism. These are low amplitude, high 40 
frequency effects, causing segment motions of 0.01-0.08 
microns peak at frequencies ranging from 0.5-500 Hz. 
They are not affected by the figure maintenance con- 
troller, being beyond its effective bandwidth. They are 
passed through directly to the panels causing acceptable 45 
levels of image jitter. 

The objective of the FMMT figure maintenance con- 
troller is to maximize the optical quality of the entire 
telescope with respect to primary mirror segment posi- 
tion errors. It does this by moving the primary mirror 50 
segments to cancel optical aberrations induced on the 
wavefront by rotational and translational disturbances 
of the segments. Our controller uses an accurate mathe- 
matical model of the FMMT optics to compute the 
control that minimizes the wavefront error of the tele- 55 
scope. This differs from previous work, w f hich seeks to 
minimize segment position errors directly, without ac- 
counting for their effect on the ultimate imaging perfor- 
mance of the telescope. 

We treat the FMMT as a quasi-static system, which is 60 
appropriate given the long time-scale of the thermal 
disturbances as compared to the maintenance control 
bandwidth. We consider only the six outer segments, 
using the center segment as a reference for the rest. We 
neglect the effects of segment surface distortions and 65 
motions of other optical elements, though these can 
easily be included in the analysis if desired. We concen- 
trate on identifying the basic limiters of performance for 
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figure control systems: controllability, observability, 
actuator and sensor noise. We address these and other 
issues in the context of the FMMT, but the methods we 
use, and the conclusions we draw, can be applied to any 
segmented-mirror telescope. 

DESCRIPTION OF THE FMMT 

Referring to FIG. 2, the FMMT main optics form an 
f/10 Cassegrain telescope, consisting of a very fast 
//0.4 segmented parabolic primary mirror 110 and a 
hyperbolic secondary mirror 120 . The FMMT also uses 
1:1 relay optics to reimage the light in a cryogenic 
dewar which contains the pixillated science detectors 
(not shown). Currently, we compute the pupil function 
of the system at the spherical reference surface (RS) 130 
located just past the secondary mirror 120 . The refer- 
ence surface 130 is aligned with the nominal spherical 
wavefront 140 of the unperturbed system. One of the 
baseline science detectors in the focal plane 142 (FP) is 
a 10X10 bolometer array which is sensitive in the 
100-300 fi,m wavelengths of light. This bolometer is 
used by the initialization controller in the initial phasing 
of the panels, as well as to collect astronomical science 
data. 

The primary mirror 110 is supported by a trusswork 
backplane structure (not shown) constructed of high 
stiffness, low thermal coefficient-of-expansion members 
and consists of mirror segments Si, S2, S3> S4, S5, S6and 
S7. These segments attach to the nodes of the truss and 
are influenced by motions of the nodes. The thermal 
and other slow motions of the optics are expected to be 
between 10-235 \xm and 10-200 /xrad over the mainte- 
nance cycle. 

Each segment panel S2, S3, S4, S5, S6 and S7 is sepa- 
rately actuated in three degrees of freedom. This is 
accomplished using three linear actuators 100 for each 
segment panel acting in the segment z-axis at the three 
segment attach points. For this analysis we resolve the 
three available control degrees of freedom into tip (ro- 
tation about the local y-axis); tilt (rotation about x); and 
piston (translation along the local z-axis) commands 
acting at the center of the panels. The segment actuator 
coordinates are centered in each panel with the x-axis 
pointing radially away and tangent to the surface, the 
z-axis along the normal, and the y-axis forming a tan- 
gent in the circumferential direction (FIGS. 3 and 4 ). 
One-sigma actuator noise levels are expected to result in 
approximately 350 nm error in piston and 350 nrad in tip 
and tilt. The three uncontrolled degrees of freedom per 
panel consist of x- and y-decenter (translation along x or 
y) and twist (rotation about z). 

Also indicated on FIG. 3 are the locations of figure 
sensors for one candidate sensor configuration. Two 
types of figure sensors are used in this analysis. Edge- 
sensors IOO0 use an array of four laser interferometers 
102 to directly measure relative piston (z-axis transla- 
tion) and hinge angle (x-axis rotation) between two 
segments, as indicated on FIG. 5. Gap sensors 1006 
using interferometric, inductive or capacitive devices 
can be used to measure relative displacements in the 
y-axis direction. One-sigma measurement noise of about 
1 nm in piston and gap was assumed; hinge angle noise 
of about 1 nrad was also assumed. 

MODELLING 

The first step in creating a controller of this type is to 
formulate the objective function in terms of the error 
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states and the control variables. In our case, the objec- 
tive function is wavefront error squared at the exit pu- 
pil, the error states are the 6 degrees of freedom (DOF) 
rotational and translational motions of the six outer 
segments, and the controls are the 3 DOF (tilt, tip and 5 
piston) actuators of each segment. The second step is to 
take the derivatives of the objective function with re- 
spect to the error states and control variables. These 
derivatives provide a linear model of the influence of 
the errors and controls on the objective. The third step 10 
is to invert the linear model (in the least-squares sense) 
to obtain expressions for the control in terms of the 
errors, as weighted by the error contribution to the 
objective function. The result is the optimal control 
gain, i.e. the control that minimizes the objective func- 15 
tion. 

A similar process is followed to create an estimator 
for processing segment figure sensor measurements into 
segment error state estimates. Here we start by express- 
ing the sensor measurements as a linear function of the 20 
segment state errors using the kinematics of the FMMT 
structure. This measurement equation is linearized and 
inverted (using weighted least-squares) to obtain equa- 
tions for estimator gains, giving segment state errors as 
a linear function of sensor measurements. 25 

In most of this analysis we treat the FMMT as a 
quasi-static system. We consider the errors to be initial- 
condition errors and we neglect structural and thermal 
dynamics. In the case of the thermal disturbances, the 
time scale is so slow compared to segment servo dynam- 30 
ics. In the case of the thermal disturbances, the time 
scale is so slow compared to segment servo dynamics 
that the controller is well approximated as a periodic 
single-stage controller without dynamics. As for the 
structural vibrational disturbances, they occur so fast 35 
compared to the maintenance controller bandwidth that 
they pass through unimpeded. 

The wavefront error is best modelled using ray-trace 
optics. Ray-tracing treats light beams as bundles of 
straight-line rays. Each ray is broken into segments, 40 
with direction changes occurring at reflective and re- 
fractive surfaces. For the FMMT, we trace a bundle of 
n (a large number) collimated rays starting at the en- 
trance pupil. This represents the light from a point ob- 
ject at infinity (not a band approximation to a distant 45 
star). We follow each ray as it hits one of the primary 
mirror segments, then the secondary mirror, and then 
finally the spherical reference surface, which has its 
center of curvature at the nominal focal point. In the 
nominal case, the pathlength of each ray is exactly the 50 
same. When the system is perturbed, the length of some 
of the rays will change. Differencing the perturbed and 
nominal pathlengths gives the ray OPDs (optical path 
differences) at the reference surface. The OPDs in turn 
determine the phase of the complex pupil function at the 55 
reference surface. The focal plane image is determined 
by propagating the wavefront determined by the pupil 
function to the focal plane using physical optics tech- 
niques. System image quality is indicated by the .RMS 
wavefront error computed directly from the OPDs. 60 

To compute the ray-trace models we follow the anal- 
ysis of Redding et al., “Linearized Ray-Trace Analy- 
sis,” International Lens Design Conference , June, 1990, 
which provides a full ray-trace analysis that can be 
performed inline with other computer programs, or 65 
separately on a stand-alone basis. We use a linearized 
model of the FMMT set up by directly calculating the 
partial derivatives of the rays at the reference surface 
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with respect to the segment states and control actions, 
to determine controllability and observability of the 
controller, and to perform covariance analyses. The 
mathematics are summarized briefly in the next few 
paragraphs. 

The analysis represents the i-1 ray segment in terms of 
its starting point pi to the vertex of the \ th surface, its 
direction r /./, and its length L/.i(see FIG. 6). Starting at 
py, the ray proceeds in the direction r/./ until it hits the 
\ th surface. Conic surfaces are represented by a dyad M 
and a vector No, which are functions of the geometric 
focal length f, eccentricity e and principal axis direction 
* 


M-(I-e 2 W) 

0) 


N 0 =-fy 

(2) 


The distance travelled to the surface is given by L /, 
which is a solution of 


+ 2 No)Lj+p]-(M- 

'pi-T 2 Aq) = 0 (3) 


The point of incidence on the surface is given by 

pJ-pl+Liru (4) 

Pi is the starting point for the next (i' A ) ray segment: the 
process repeats until the ray reaches the last surface. 
The direction of the ray reflected from the surface is 
given by 


r=R-i 


(5) 


where R is the “reflection dyadic” 

R=l— 2NN 


( 6 ) 


and the normal at the point of incidence is 


N=unit (No+Mpb (7) 

The effect of perturbations either in the incoming ray 
or in the surface attitude and position is to change the 
outgoing ray state. Perturbed rays are represented in 
terms of small ray direction perturbations dr, transverse 
aberration or “beam walk” y, and the optical path differ- 
ence (OPD) dL (FIG. 7). The Redding et al. publication 
presents closed-form equations for the derivatives of 
ray states with respect to perturbations of the incident 
ray and of the surface perturbations 6 (rotation) and 8 
(translation). These derivatives are chained to develop a 
linearized model of the complete optical system. For 
the FMMT we are concerned with the OPD at the exit 
pupil only. The linear model takes the form 

wo = Cxo (8) 

Here e,rar/w/ o is a vector of n ray OPDs, and x is the 
error state, composed of six 6-v zctor $_x S egi (one for each 
of the outer segments). We take the x S egi s to be in seg- 
ment coordinates, so that each is in the form 

x s££/ = [ $ z] ^ • W 

B x , By and B z are rotations about the x, y and z axes 
while 8*, d y and 8 2 are rotations about the x, y and z axes 
of the corresponding segment. 
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C is thus an array of partial derivatives relating the 
ray OPDs to element perturbations: 


the two segments: relative piston (local z translation) 
and relative tilt (local x or y rotation), as sketched in 
FIG. 5. The gap sensors 1006 are capacitive devices 
that, when used, are mounted in pairs at the same loca- 
tions as the edge sensors. They provide measurements 
of relative gap (local x or y translation) (FIG. 5). Both 
sensor systems can provide excellent performance, 
yielding accuracies in the nanometers and nanoradians. 

Direct derivation of the linearized measurement 
equations is straight-forward. In general, the form of the 
measurement equation is 


Two other linear models are required in the deriva- 
tion of the controller. The first is the control action 
matrix A, which scales actuator motions u directly into 15 
segment coordinates. A is a 36 by 19 matrix which 
simply projects the three controls (tig 6 X > tilt $ y and 
piston 8 Z ) into segment coordinates, u is the control 
vector made up of the six individual segment controller 
vectors u^/, where each u se gi is 20 


'u'seqi = °y 


Given that the controls act along the segment coordi- 
nates, the 6 by 3 A matrix for one segment is 


x segi Usegi 0 
0 

' 0 0 
0 0 Zsegi 


As an example, consider the edge and gap sensors be- 
tween two panels i and j as sketched on FIG. 8. The 3 
measurements are formed into the measurement vector 
z ^ which is a function of the two segment states 


hinge angle 
piston 
gap 


bzy bzjj Iran 

»*j JU J 


The full A matrix is 


25 

The partials are 

derived by inspection as 



-Xjj 

0 0 0 

(19) 

30 

3*1/ 

-zT- 

z u 



Pip? 

y v 



and 



35 

-*9 

0 0 0 

(20) 


13T = 

z y 



yW 

y ‘J 



A S eg2 0 0 


The full set of measurement equations for each of the 
configurations is an elaboration of this simple case, with 
the measurement vector z consisting of the z# from each 
sensor; and xq being the full error state 


0 0 A se g-j 

The effect of actuator motions in segment state coordi- 
nates is 


z=Hzo + r 

Here r is sensor noise. The measurement matrix H is 
built up of the partials of the measurements as above 


x oc/ =Au (14) 

Adding the effect of actuator motions to those of the 
segment errors, the total wavefront error is 

1*1 — wq+CAu-Cxq~\- CAu (15) 

The controller objective function is the total wavefront 
error squared: 


**12 **12 
dx\ 35$ 


3*67 3*67 

5*6 3*7 


J — w J j 


(16) 60 WAVEFRONT CONTROLLABILITY 


The last part of the modeling task is the derivation of 
the measurement equations. These depend on the partic- 
ular sensor configuration, and we will consider several. 

All use the same basic elements: edge sensors 100a and 65 
gap sensors 1006 located between adjacent panels (FIG. 

5). The edge sensors 100a use laser interferometers 
configured to measure two degrees of freedom between 


In keeping with our use of quasi-static model of the 
FMMT, we treat the control problem as a single-stage 
optimal control problem. We derive a control that takes 
the objective function— the wavefront error 
squared — to its minimum possible value in a single in- 
stantaneous step. This is a good approximation for our 
system, as the time scale of the control is very short 
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compared to the thermal disturbances. In the real sys- 
tem the control will be implemented using a dynamic 
controller that takes the panel and actuator dynamics 
into account. Our intention here is to address the under- 
lying controllability question: Can the segments be con- 5 
trolled to completely eliminate the effects of segments 
error on the wavefront? The answer is a function of the 
optics and the system geometry, and is the same for 
both the quasi-static and dynamic controllers. 

The control problem is to minimize _the objective jq 
function J over the transition from wo to wj. We derive 
the optimal control by finding the value of u that ren- 
ders J stationary, or (in other words) satisfies dJ =0 for 
non-zero values of du: 

dJ=dw\^\=dpA T C T {CAi?+Cx q )^0 (22) 15 

This condition is satisfied by the control 

uJ p ,= _[A r C r CA]“ 1 A 7 'C 7 '05 (23) 

20 

The control as implemented includes actuator noise q, 
and uses the estimated panel error state !«,. Defining G 
as the control gain matrix, the control as implemented is 

~?=-Gxest+7 (24) 25 

where 

G = [A T C T CA]~ l A T C T C (25) 

The same result can be derived by choosing G to mini- 30 

mize the variance of the residual wavefront error, or as 
a minimization with the transition equation 24 adjoined 
to the objective function J as an equality constraint. 

Consider the general problem of estimating the seg- 
ment error^states. We start with some initial segment 35 
error state xo_with covariance Xo- The sensor measure- 
ment signals z are obtained via the measurement equa- 
tion discussed above 

T=m 8+7* (26) ^ 

where r is the measurement noise, characterized as 
zero-mean noise with covariance R. The particulars of 
H depend on which sensor configuration is used, as 
discussed previously. We assume that the previous best 
estimate of the segment state was zero. The objective 45 
function for the weighted least-squares estimate of xo is 

J=a){7 T X 0 ~^+{7-H7o) T R-H^H7o))( (27) 

By minimizing the wavefront error, the controller 50 
maximizes the overall optical quality of the telescope 
with respect to the modelled degrees of freedom. Other 
error sources, such as motions of the secondary mirror, 
relay optics and focal plane, and thermal deformations 
of the primary mirror segments, also contribute to the 55 
overall system performance. The framework can easily 
be expanded to include all these effects, provided ap- 
propriate sensors and actuators are added as well. To 
realize this, the state vector, which in this paper, con- 
tained only the six rigid body states of each of six pri- 60 
mary reflector panels, would be augmented to include 
the six secondary mirror rigid body states, the six focal 
plane rigid body states, and as many primary segment 
deformation states (e.g., thermal mode shapes) as are 
deemed important. Similarly, the control vector would 65 
be augmented with secondary mirror tip, tilt, x-decen- 
ter, y-decenter and focus, and panel deformation con- 
trol inputs. The measurement vector would be modified 
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to include more samples of the primary segment sur- 
faces, secondary positions and orientations, and focal 
plane positions and orientations. The basic mathematics 
of the problem, which yield an optimal control solution 
for the modelled degrees of freedom, are unchanged 
from the results presented in this paper. 

The controller should outperform any prior art ap- 
proach with respect to wavefront error. Telescope opti- 
cal parameters, such as focal ratio and relative segment 
size, will significantly affect the performance of any 
controller. 

While the invention has been described in detail with 
specific reference to preferred embodiments thereof, it 
is understood that variations and modifications may be 
made without departing from the true scope and spirit 
of the invention. 

What is claimed is: 

1. In a complex optical system having a light-receiv- 
ing aperture, an objective pupil and plural optical sur- 
faces of respective shapes disposed in certain locations 
and orientations which determine an optical path for 
each ray in a bundle of rays travelling from said aper- 
ture to said objective pupil, said shapes, locations and 
orientations being subject to perturbations due to exter- 
nal influences from ideal shapes, locations and orienta- 
tions, respectively, said perturbation corresponding to 
changes in the optical path of each of said rays, said 
system further having plural actuator means associated 
with said plural optical surfaces for changing said loca- 
tions and orientations in response to an actuator com- 
mand vector specifying plural commands for said plural 
actuator means and sensor means for sensing each com- 
ponent of perturbation of the location and orientation of 
each of said plural optical surfaces, a feedback control 
system for compensating for said perturbations, said 
feedback control system comprising: 

means for generating from said ideal shapes, locations 
and orientations a sensitivity matrix which defines, 
for each ray in said bundle of rays, a sensitivity of 
its optical path to each component of perturbation 
of the location and orientation of each of said plural 
optical surfaces, said sensitivity representing the 
rate of change of said optical path responsive to a 
corresponding component of perturbation; 
means for generating an actuator response matrix 
which relates said actuator command vector to 
corresponding perturbations of said locations and 
orientations caused by said actuator means in re- 
sponse to said actuator command vector; 
means for generating from said sensitivity matrix and 
from said actuator response matrix a gain matrix by 
determining a condition whereby a combination of 
said changes in the optical paths of all of said rays 
is minimized at the objective pupil; and 
means for producing an optimum actuator control 
vector by multiplying by said gain matrix an error 
vector comprising each component of perturbation 
in said shapes, locations and orientations of said 
plural optical surfaces sensed by said means for 
sensing, and transmitting commands corresponding 
to said optimum actuator control vector to respec- 
tive ones of said actuator means so as to control the 
shapes, locations and orientations of said plural 
optical surfaces. 

2. The system of claim 1 further comprising: 
means connected to said means for sensing for pro- 
viding modifications to said ideal locations and 
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orientations attributable to perturbations in said 
locations and orientations sensed after said plural 
actuator means have responded to said actuator 
control vector; and 

wherein said means for generating comprises means 5 
for modifying said sensitivity matrix in response to 
said modifications to said locations and orienta- 
tions. 

3. The system of claim 1 wherein said perturbations 
are defined with respect to orthogonal basis vectors and 10 
wherein said sensitivity matrix is an orthogonal matrix 
C, said actuator response matrixes an orthogonal matrix 
A, said error vector is a vector x, said actuator control 
vector is a vector u, and wherein said means for gener- 
ating a gain matrix comprises means for determining 
how to minimize a total wavefront error 
[C(x 4- Au] ^CXx -f Au] . 

4. The system of claim 3 wherein said means for pro- 
ducing a gain matrix produces the following gain ma- 
trix: [ATCTCA]-l[ATCTC]. 

5. A complex optical system having a light-receiving 
aperture, an objective pupil and plural optical surfaces 
of respective shapes disposed in certain locations and 
orientations which determine an optical path for each 
ray in a bundle of rays travelling from said aperture to 
said objective pupil, said shapes, locations and orienta- 
tions being subject to perturbations due to external 
influences from ideal shapes, locations and orientations, 
respectively, said perturbation corresponding to 
changes in the optical path of each of said rays, said 
system comprising: 

plural actuator means associated with said plural 
optical surfaces for changing said locations and 
orientations in response to an actuator comand 
vector specifying plural commands for said plural 
actuator means; 

means for sensing each component of perturbation of 
the location and orientation of each of said plural 
optical surfaces; 

means for generating from said ideal shapes, locations 40 
and orientations a sensitivity matrix which defines, 
for each ray in said bundle of rays, a sensitivity of 
its optical path to each component of perturbation 
of the location and orientation of each of said plural 
optical surfaces, said sensitivity representing the 45 
rate of change of said optical path responsive to a 
corresponding component of perturbation; 

means for generating an actuator response matrix 
which relates said actuator command vector to 
corresponding perturbations of said locations and 50 
orientations caused by said acutator means in re- 
sponse to said actuator command vector; 

means for generating from said sensitivity matrix and 
from said actuator response matrix a gain matrix by 
determining a condition whereby a combination of 55 
said changes in the optical paths of all of said rays 
is minimized at the objective pupil; and 

means for producing an optimum actuator control 
vector by multiplying by said gain matrix an error 
vector comprising each component of perturbation 60 
in said shapes, locations and orientations of said 
plural optical surfaces sensed by said means for 
sensing, and transmitting commands corresponding 
to said optimum actuator control vector to respec- 
tive ones of said actuator means so as to control the 65 
shapes, locations and orientations of said plural 
optical surfaces. 

6. The system of claim 5 further comprising: 
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means connected to said means for sensing for pro- 
viding modifications to said ideal locations and 
orientations attributable to perturbations in said 
locations and orientations sensed after said plural 
actuator means have responded to said actuator 
control vector; and 

wherein said means for generating comprises means 
for modifying said sensitivity matrix in response to 
said modifications to said locations and orienta- 
tions. 

7. The system of claim 5 wherein said perturbations 
are defined with respect to orthogonal basis vectors and 
wherein said sensitivity matrix is an orthogonal matrix 
C, said actuator response matrixes an orthogonal matrix 
15 A, said error vector is a vector x, said actuator control 
vector is a vector u, and wherein said means for gener- 
ating a gain matrix comprises means for determining 
how to minimize a total wavefront error 
[C(x -f Au] 7lC(i + Au] , 

20 8. The system of claim 7 wherein said means for pro- 

ducing a gain matrix produces the following gain ma- 
trix: [A A] - 1 [A r C ^C] . 

9. A method of maintaining a complex optical system 
having a light-receiving aperture, an objective pupil and 
25 plural optical surfaces of respective shapes disposed in 
certain locations and orientations which determine an 
optical path for each ray in a bundle of rays travelling 
from said aperture to said objective pupil, said shapes, 
locations and orientations being subject to perturbations 
30 due to external influences from ideal shapes, locations 
and orientations, respectively, said perturbation corre- 
sponding to changes in the optical path of each of said 
rays, said system having plural actuator means associ- 
ated with said plural optical surfaces for changing said 
35 locations and orientations in response to an actuator 
command vector specifying plural commands for said 
plural actuator means: 

generating from said ideal shapes, locations and ori- 
entations a sensitivity matrix which defines, for 
each ray in said bundle of rays, a sensitivity of its 
optical path to each component of perturbation of 
the location and orientation of each of said plural 
optical surfaces, said sensitivity representing the 
rate .of change of said optical path responsive to a 
corresponding component of perturbation; 
generating an actuator response matrix which relates 
said actuator command vector to corresponding 
perturbations of said locations and orientations 
caused by said actuator means in response to said 
actuator command vector; 
generating from said sensitivity matrix and from said 
actuator response matrix a gain matrix by deter- 
mining a condition whereby a combination of said 
changes in the optical paths of all of said rays is 
minimized at the objective pupil; 
sensing each component of perturbation of the loca- 
tion and orientation of each of said plural optical 
surfaces; 

producing an optimum actuator control vector by 
multiplying by said gain matrix an error vector 
comprising each component of perturbation in said 
shapes, locations and orientations of said plural 
optical surfaces sensed by said means for sensing, 
and transmitting commands corresponding to said 
optimum actuator control vector to respective ones 
of said actuator means so as to control the shapes, 
locations and orientations of said plural optical 
surfaces. 
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10 . The method of claim 9 further comprising: 
computing, from perturbations in said locations and 

orientations sensed after said plural actuator means 
have responded to said actuator control vector, 
modifications to said ideal locations and orienta- 5 
tions; and 

modifying said sensitivity matrix in response to said 
modifications to said locations and orientations. 

11 . The method of claim 9 wherein said perturbations 
are defined with respect to orthogonal basis vectors and 10 
wherein said sensitivity matrix is an orthogonal matrix 
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C, said actuator response matrixes an orthogonal matrix 
A, said error vector is a vector x, said actuator control 
vector is a vector ii, and wherein said step of generating 
a gain matrix comprises determining how to minimize a 
total wavefront error [C(x+ Aii] 7 lC(x + Au]. 

12 . The method of claim 11 wherein said step of pro- 
ducing a gain matrix produces the following gain ma- 
trix: 

[A A] ~ [ A ^C^C]. 

* * * * * 


15 


20 


25 


30 


35 


40 


45 


50 


55 


60 


65 



